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With exposure to air, many liquid metals spontaneously generate an oxide layer on their surface. 
In oscillatory rheological tests, this skin is found to introduce a yield stress that typically dominates 
the elastic response but can be tuned by exposing the metal to hydrochloric acid solutions of different 
concentration. We systematically studied the normal impact of eutectic gallium-indium (eGaIn) 
drops under different oxidation conditions and show how this leads to two different dynamical 
regimes. At low impact velocity (or low Weber number), eGaIn droplets display strong recoil and 
rebound from the impacted surface when the oxide layer is removed. In addition, the degree of 
drop deformation or spreading during the impact is controlled by the oxide skin. We show that 
the scaling law known from ordinary liquids for the maximum spreading radius as a function of 
impact velocity can still be applied to the case of oxidized eGaIn if an effective Weber number 
We* is employed that uses an effective surface tension factoring in the yield stress. In contrast, 
no influence on spreading from different oxidations conditions is observed for high impact velocity. 
This suggests that the initial kinetic energy is mostly damped by bulk viscous dissipation. Results 
from both regimes can be collapsed in an impact phase diagram controlled by two variables, the 
maximum spreading factor P,n = Ro/Rm, given by the ratio of initial to maximum drop radius, and 
the impact number K = We* /Re'*'^^, which scales with the effective Weber number We* as well 
as the Reynolds number Re. The data exhibit a transition from capillary to viscous behavior at a 
critical impact number Kc ~ 0.1. 



I. INTRODUCTION 

Normal impact of liquid metals is important for a wide variety of industrial processes and applications, including 
electronic fabrication!^! , thin film coatingPl'I^l and heat conductor production W. Typically, precise control of the 
drop spreading is desired, in particular for processes that require reproducible and predictable behavior to build up 
layers from successive drop impact cvcntsl^l^l^l . However, except for mercury most liquid metals develop an oxide 
skin when exposed to air^^l'I^l, resulting in non-Newtonian behavior. One of the direct consequences is that the 
oxide layer consumes a portion of kinetic energy and thus deviates their spreading dynamics from ordinary fluids 
in the impact. Also, oxidized liquid metals can generally form non-spherical drops since the surface skin not only 
prevents contact between air and the bulk of the material but also the minimization of the surface energy. Previous 
investigations observed the effect of the oxide skin in various measurements, including rheology, wetting capability 
and surface tension I"*! ' ~ . However, a detailed experimental study of how the presence of an oxide skin controls 
the dynamics of liquid metals impinging onto a solid surface has been lacking. This is the focus of the present paper. 

In the absence of splashing[^°l~[^^l, previous experiments on Newtonian fluids mainly focused on the geometric 
deformation of the drop, such as the maximum lateral spreading distance during the impact [^^!~[^^1. Conventionally, 
this deformation is modeled from an energy conservation point of view, considering the different energy scales affecting 
the impact and spreading processes: viscous drag dissipates the momentum inside the liquid metal: surface tension 
stores the initial kinetic energy; and friction from the substrate resists the spreading. We can then associate different 
regimes with conditions under which certain energy scales dominate the dynamics. For instance, Clanet et al. I^^l 
showed experimentally for Newtonian fluids that the maximum spreading radius scales as a power of the Weber 
number in the capillary regime and a (different) power of the Reynolds number in the viscous regime. Here, we will 
show that similar scaling laws can still be used for liquid metals if the relation between skin stress and oxidation 
condition is appropriately taken care of. 



To realize a systematic control of the skin effect, we use eutectic gallium-indium (eGaIn) as a model fluid, for which 
oxidation can be managed by immersing the metal in a bath of acid. The acid not only prevents continued oxidation 
by mitigating the contact with air, but also initiates a chemical reduction process that, in equilibrium, produces 
a skin thickness corresponding to the particular reaction conditions. Therefore, we can tune the skin by varying 
the acid concentration and use this to quantitatively analyze the role of the oxide layer in the impact dynamics. 
The purpose of the paper is to lay down an experimental framework for describing the impact behavior of liquid 
metals. The remainder of the paper is organized as follows. We start with the rheological characterization (Sec. 
Ill) to confirm the change of skin-induced stresses under different acid concentrations. Then we qualitatively show 
how the change of surface elasticity affects the impact behavior on glass (Sec. IV). Within the capillary regime, we 
demonstrate how unoxidized eGaIn drops can rebound from the substrate at very low Weber numbers (Sec. Va). 
In order to collapse the spreading data for various surface; conditions, we introduce an effective Weber number that 
depends on drop size and skin stress (Sec. Vb). We then show that the skin effect does not play a significant role 
in the high- velocity, viscosity-dominated regime (Sec. VI). Finally, we develop a phase diagram that collapses data 
from both regimes by using appropriately scaled dimensionless groups (Sec. VII), before concluding in Section VIII. 

II. EXPERIMENTAL 
A. Materials 

We purchased the eutectic gallium-indium (eGaIn, purity 99.995%) from Sigma-Aldrich. This material has a 
melting point around 15°C. All experiments were performed at room temperature (~ 24^C ), thus the materials 
stayed in its liquid state at all times. At this temperature, liquid cGaln has a density p — 6.25g/ml and a nominal 
viscosity rji = 1.92 x 10~^Pa- s. For the acid bath, we use hydrochloric acid (HCl) of varying concentration. As 
substrate for the impact experiments, we used a corrosion-resistant glass plate of 8mm thick. 

B. Rheological Measurements 

To characterize the stress response of eGaIn, we performed rheological test in an Anton Paar MCR 301 rheometer 
with parallel-plate geometry. Figure 1 illustrates our experimental setup. In the measurement, the top rotating tool 
had a radius r of either 12.7mm or 4.1mm. The cGaln sampls were settled in the gap between rotating top and fixed 
bottom plate. The gap size d was varied from 1 to 2mm. 

For unidirectional shear measurements, the rheometer shears the fluid by rotating the tool at a constant angular 
rotation rate u, and the torque T required to shear the fluid is measured in the steady state. The average shear 
stress is calculated as r = 2T/7rr^ and the average shear rate at the plate edge is 7 = ur/d. 

For oscillatory tests, on the other hand, the rheometer tool oscillates back and forth at a fixed frequency and strain 
amplitude 7. The storage modulus, defined as G' = r/7, provides an indicator of the elastic response to the shear. 

To contain both acid bath and sample during the measurements, a cup was constructed consisting of a clear acrylic 
cylinder with 31.75mm inside diameter and glued to a titanium plate of thickness 0.51mm that formed the bottom. 
Further, the transparency of the cup provides the opportunity to monitor the sample underneath the tool. Due to 
the high surface energy with most solid surfaces, eGaIn does not stick to the titanium plate. For this reason, we fill 
the eGaIn sample to a height larger than the gap size d so there is no extra space for the liquid to spill outward. 
Also, a 13mm thick acrylic cylinder was glued to the bottom of the top plate of the rheometer to protect it from the 
acid (see Fig. 1(a)). From measurements at low shear rate, the stress resolution of our rheometer was found to be 
0.06Pa. All rheological measurements were performed multiple times to check for repeatability. Because the standard 
deviation from those repeated tests was 5% or better, in the following we show data from only one measurement for 
clarity. 
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C. Imaging system 



The imaging setup for the drop impact is sketched in Fig. lb. Drops were slowly extruded from a syringe with steel 
nozzle (inner radius of 0.48mm and 4.1mm) by using the RAZEL 99E syringe pump (the pump rate was fixed at 
20ml/hr). The fahing height was varied from 1 to 200cm to give an impact velocity Vo — (0.4 ^ 6.3) ±0.15m/s to the 
drop before impinging on the glass substrate. Using an inclined reflective mirror under the substrate, we were able 
to monitor the entire process from both side and bottom. To highlight the drop profile, the drops were backlit by 
two white light sources (12V/200W, Dedolight). We used a Vision Research Phantom vl2 camera with a macro lens 
(Nikon Micro 105mm) to video-capture the impact process. The frame rate reached 17, OOOfps at 380 x 540 pixels 
with a spatial resolution around 15/im. The camera was adjusted to include both the side and bottom images (from 
the mirror) in the same field of view. Fig. Ic shows typical images captured by the camera. 
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FIG. 1: (a). Sketch of the rheometer setup. An acrylic cup (31.75mm inner diameter and 13.10mm height) was constructed 
to contain the acid bath, in which the rotating tool shears the eGaIn sample, (b). Schematic representation of our imaging 
setup. The syringe pump can vertically move along the trail by two meters that is made of four steel rods. Drops during the 
impact can be lit from the back or the top, depending on the perspective of the video taken by a Phantom vl2 camera. The 
highest frame rate is up to ~ 10^ fps. (c). Two typical images of a liquid eGaIn drop impacting a smooth glass substrate, 
showing side and bottom views simultaneously captured with the help of the mirror shown in panel (b). 

For a given nozzle diameter, the initial drop radius was found to be reproducible with 5% uncertainty or less. To 
quantify the relative deformation of the drops with time t, we use the spreading factor P, defined as P = R{t)/Ro, 
where R{t) is the spreading radius obtained by averaging the distance from the advancing contact line to the center 
of the initial impact and the initial drop radius Rq is measured from the half-width of the horizontal cross-diameter 
of the drop before impact. All the droplet-related geometric dimensions reported in the following were extracted 
from the videos using edge-detection algorithms in Matlab. 



D. Sample Preparation 



Prior to all rheological measurements, the samples were pre-sheared, or washed, at 7 = 60s ^ for 10 mins while 
immersed in the HCL solution, so that the skin sufhcicntly reacted with the acid. The eventual chemical equilibrium 
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determines the amount of oxide left on the surface of eGaln. The steady state was confirmed by slowly ramping the 
shear rate first up and then down, and checking there was no hysteresis in the measured stress, except when eGaln 
was exposed to air, in which case oxidation is not reversible under shear^^l . 

Samples used in our study of the impact process were prepared in the same way and also sheared at constant rate of 
60s^^ for 10 minutes. After this washing step, we used a plastic syringe with a steel nozzle tip to extract eGaln from 
the bath. Then the syringe was inserted into the pump and liquid eGaln was slowly ejected from the nozzle to form 
the drops. Since the volume of a single droplet was around 0.15ml, it took approximately 0.15ml/(20ml/hr) = 20s 
to fully extrude the drop until it fell off. Together with the time spent in mounting the syringe, impact occurred in 
less than 30s after the washing step. Also, to avoid possible oxidation during the time of performing multiple impact 
experiments, the data we used is usually from the first drop extruded from the nozzle. Therefore, fresh oxidation, 
which usually initiates after a sample has been exposed to air for a couple of minutes, can be neglected for our 
measurements . 



III. RHEOLOGY RESULTS 
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FIG. 2: Shear stress vs. shear rate curves for liquid eGaln in HCl baths with different acid concentrations. A detailed discussion 
of the similar result can be referred to [7] (Sec. III. A. 1. , Fig. 4). 

The drop impact dynamics of non-Newtonian fluids are usually modeled based on their intrinsic rheological prop- 
erties. However, we have shown previously that surface oxidation is crucial for the shear stress measurement of eGaln 
in the steady state'^l : for instance. Fig. 2 provides plots of shear stress r against shear rate 7 for eGaln submerged 
in acid baths of different concentrations. Initially, 2ml of eGaln was directly placed on the bottom plate of the 
rheometer and exposed to air. Such an oxidized droplet appears dirty and wrinkled on the surface, while the bulk 
appears shiny and mirror-like if we slice open the skin. When exposed to air (red upward pointing triangle points 
in Fig. 2), eGaln displays a significant yield stress Ty ^ lO^Pa, indicating an effective solidification of the material. 
Adding HCl into the bath can eliminate the skin effect. As shown in Fig. 2, when the acid concentration reaches 
Chci > O.IM, Ty dramatically drops by four orders of magnitude to nearly zero within the rheometer resolution. 

This vanishing yield stress at sufficiently high acid concentration suggests that pure eGaln behaves as an ordinary 
viscous fluid, for which the dynamic viscosity is given by the ratio between shear stress and shear rate in the laminar 
regime. However, because of the large density of eGaln, inertial effects can easily play a more essential role in the 
shear flow than they would for normal fluids. As a result, at high shear rate all traces in Fig. 2 collapse onto the 
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typical scaling of incrtial flow[''l'[^°l, r ~ 7^/^ (solid line). 

Fig. 2 confirms that the measured yield stress of the eGaIn samples in air is associated with oxidation. The relative 
contributions to the overall stress from the bulk of the material and from the skin can be extracted from measurements 
on samples of different sizes. To this end, we did oscillatory viscoelasticity measurements using two different sizes 
for the rotation plate, 8.2mm and 25.4mm in diameter, with correspondingly different sample volumes and exposed 
surfaces around the perimeter of the plate. 

Performing strain oscillations at fixed angular frequency uj = 0.5 rad/s and gap size d = 1mm, the applied average 
strain 7 was ramped from 0.001 to 10. For the oxidized eGaIn sample, over the range of 7 < 0.1, the elastic modulus 
G' of the oxidized eGaIn was found to be two orders of magnitude larger than the loss modulus G" (this is similar 
to the observation by Larsen et al.I^l, except that they used uj = 1 rad/s). 

For ordinary viscoelastic materials, whose stress originates from the elastic response of the bulk, G' is supposed 
to be an intrinsic quantity, independent of sample size. However, the apparent elastic modulus of oxidized eGaIn is 
found to be more complicated. Fig. 3a shows G' against strain 7 at different acid concentrations. When the sample 
is exposed to air (red circles) or O.OIM acid solution (blue circles), the plot shows significant size dependence: G' is 
reduced by a factor around three when the plate diameter is increased by the same factor. In addition, a decrease 
of G" starts to occur at 7 ~ 10~^, which corresponds to the end of the linear elastic regime and signals breaking 
of the skin. By contrast, for Chci > O.IM (pink and green circles), the obvious size dependence of G' disappears. 
Instead, over the range of 7 scanned, G' constantly stays below the rheometer resolution limit (~ 0.06Pa, the dashed 
black line in the figure), which indicates a very small bulk modulus of liquid eGaln. Thus, the behavior of G' more 
consistent with ordinary fluids is recovered at high acid concentration, or weak skin effect. In other words, G' may 
not be an appropriate, intrinsic parameter when the skin exists. 

We can model the sample size dependence by considering the elastic response of the sample as originating from 
two sources, the bulk and the skin, and treating the skin as a very thin, effectively two-dimensional layer. Inside 
such layer, in direct analogy to a bulk stress, the surface stress relates the in-plane surface strain 7s to the surface 
elastic modulus G'^ by 

Ts = • 7. (1) 

Here, Tg and G'^ are defined as forces per unit length, i.e., represent a surface tension, and are intrinsic quantities 
characterizing the surface. Generally, the total measured torque T incorporates contributions from both the skin 
and the bulk of the material, so that 

T = 2nr\G',j) + Inr^iG',^) (2) 

GJi is the regular bulk modulus. Also, we assume that there is no relative displacement between the skin and bulk 
so that 7s = 7. Using t = 2T/nr^ and defining the measured, effective elastic modulus of the sample as G' = T/7, 
we have 

G'=4G'Jr + G', (3) 

Therefore, when the surface elasticity dominates the shear response [G'^/r ^ G[,), the measured modulus G' can be 
directly related to the surface modulus by G'^ = [G' ■ r)/4.I^l 

Fig. 3b shows the same data as Fig. 3a. but rciscaled by multiplying with (r/4). Now the traces for different plate 
diameters collapse at low acid concentrations (zero or O.OIM HCl), while they are still kept below the resolution 
(dashed line) at high acid concentrations {Chci)- This indicates that indeed the surface oxidation dominates at 
low acid concentration and that the observed size dependence of the measured modulus G' in Fig. 3a reflects the 
properties of the skin. In addition, since the data is reversible in different acid concentrations, the skin effect is 
localized on the surface and does not affect the properties of the bulk. Together, Figs. 3a and b demonstrate that 
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• NoHCI(25mm) • 0.01 M {25mm) • 0.1 M (25mm) . 0.2M (25mm) 
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FIG. 3: Apparent bulk (a) and surface (b) elastic moduli as a function of applied average strain for eGaIn under uj = 0.5Hz 
oscillatory tests. Different colors (red, blue, pink, green) indicate experiments performed at various acid concentrations. Plates 
with two diameters, 25mm (solid circle) and 8mm (open circle), were used in the measurements. The dotted straight lines 
represent the rheometer resolution. 

the tiny bulk stress GJ, only shows up when the oxidation is removed. 

Therefore, for sufficiently oxidized eGaIn G'^ should be used as intrinsic parameter, while in more concentrated 
acid baths the measured G" is the appropriate intrinsic parameter, reflecting the bulk modulus of pure eGaln. 

IV. DROP IMPACT UNDER DIFFERENT OXIDATION 

By ejecting eGaln from the same nozzle and the same falling height, we generated droplets with reproducible 
impact velocity and radius, Vq — (1.02 ± 0.12)m/s and Rq — (6.25 ± 0.10)mm. Fig. 4 shows three image sequences 
for eGafn drops with different skin strengths. The left column shows the impact of an air-oxidized eGaln drop not 
pre- washed in acid. A long tail at the top end of the drop is formed when the fluid detaches from the nozzle. Different 
from ordinary liquids, the oxide skin prevents the fluid to freely relax the surface energy, so that this non-spherical 
geometry is kept unchanged during the falling stage. After the impact occurred, a thin liquid metal sheet or lamella 
expands rapidly along the smooth substrate. 

At i = 1.2ms, the lower end of the tail structure reaches the thin liquid sheet, which stimulates a secondary impact 
at the liquid-liquid interface. The induced surface wave (Fig. 4, the third image in the left column) propagates along 
the liquid sheet with the speed of shallow water waves'^^l and eventually catches up with the liquid-glass contact line 
(at 6.2ms). The surface wave is generally found in all impacts of oxidized eGaln because of the non-spherical shapes 
of the initial drop. 

In contrast, cleaning the samples with acid removes the oxide and weakens the skin effect. The middle and right 
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FIG. 4; Typical image sequences of eGaIn drops impacting onto a glass substrate. Drops are initially washed in HCl solution 
as indicated. For all three images sequences shown above, the impact velocity was kept at Vq = (1.02 ± 0.12)m/s and the 
initial drop diameter was Ro = (6.25 ± 0.10)mm. 



columns in Fig. 4 show images of drops pre-washed in O.OIM and 0.2M HCl respectively. Since the yield stress is 
reduced to ^ lOPa at O.OIM HCl, the length of the tail becomes much shorter (middle column) and only a very 
week surface wave appears {t — 3.1ms). If we keep increasing the HCl concentration to 0.2M, the acid becomes 
strong enough to fully eliminate any observable skin effect, and eGaIn shows no difference in the spreading behavior 
with ordinary liquids (right column). Generally, we do not observe any splashing of eGaIn since the surface tension 
(> 400mN/m) is much larger than ordinary liquids. 

Another major feature exhibited in Fig. 4 for a relatively large drop {Rq ^ 6mm) is the significant variance in the 
maximum spreading radius with and without oxide skin. Simple inspection of the images gives a 20% difference in the 
final stage radius between left and right column, a result that is repeatable within < 5% experimental uncertainty. 
Therefore, the skin not only affects the drop shape but also resists the spreading of the contact line along the glass 
substrate. Since the only controlling parameter is the acid concentration, the surface elasticity must play a critical 
role in determining the maximum spreading distance. 



V. CAPILLARY REGIME 

For any drop impacting on a hard surface, its kinetic energy immediately before impact will play a key role in 
the subsequent outcome. When the impact speed is low and the surface is relatively smooth, the drop typically 
spreads out and comes to rest as a truncated sphere or thin, disc-like film, depending on the equilibrium contact 
anglel^^l. When the impact velocity is sufficiently small to neglect viscous dissipation during droplet deformation, the 
balance between kinetic and surface (or capillary) energies becomes the controlling factor. In this capillary regime, 
the dimensionless Weber number, We, gives the ratio of inertial to surface stresses. Conventionally, for a drop of 
radius Rq, surface tension of a , density of p and impact velocity of Vq , the Weber number is defined as 

We = (4) 

In our experiments, due to the dramatically high surface tension of liquid eGaIn, the capillary regime is easy to 
achieve and commonly observed even at moderately large impact velocities. 



7 



No Wash 0.2M 




t (ms) 



FIG. 5: (a). Relative spreading radius R{t)/Ro as a function of time after impact for both unwashed (red) and HCl- washed 
(blue) eGaIn drops impacting a glass surface at velocity Vo = (0.83 ± 0.13)m/s. The experiments were repeated five times; the 
error bars correspond to the standard deviation (b). Image sequences for comparing the impact of a skin-covered drop (left), 
which sticks to the glass, with that of a non-oxidized drop (right), which rebounds from the surface. 

A. Rebound Control at Low Weber Number 

One of the striking outcomes that can be observed in the capillary regime is drop rebound. This is usually achieved 
by coating the substrate with a thin layer of super-hydrophobic material (such as wax or polymer)!^^!'!^-^!, in order 
to form a large contact angle (> 160*^ ) between drop and substrate. In this situation, kinetic energy is rapidly 
converted into surface energy at the drop/air interface during the spreading. After the spreading front reached its 
maximum radius, the built-up surface energy can be released by pulling the contact line backward and potentially 
lifting the drop off the surface. 

Previous studies have shown that the skin effect determines the degree to which eGaIn wets glass'^''"!. Generally, 
oxidized eGaIn wets most solid surfaces well, while pure, unoxidized eGaIn becomes perfectly non- wetting (~ 180'^). 
Intermediate contact angles can be set by the acid concentration used to wash the sample. Therefore, a precise 
rebound control of eGaIn drops can be carried out by adjusting the acid bath. 

Fig. 5 compares the behavior of drops of the same size but using oxidized and pure (washed in 0.2M HCl) eGaIn 
as they impact a glass surface with equal impact velocity Vq = (0.83 ± 0.13)m/s. At each time step, the spreading 
radius R{t) was measured and R{t)/Ro is plotted against time (Fig. 5a). For oxidized eGaIn (red upward pointing 
triangles), the drop expands rapidly along the surface until the spreading front arrives at the maximum i?,„ ^ 1.8i?o 
at t 2ms. Subsequently, even though the fluid slightly varies its surface shape to relax the interfacial energy, the 
contact line stops moving forward. Therefore, R{t) keeps its value at ^ 1.8i?o- 

The behavior of a drop of eGaIn washed in 0.2M HCl at least initially {t < 1ms) appears similar to that of the 
skin-covered drops, since spreading during this early stage is mostly due to geometric deformation. The maximum 
spreading radius, i?„i ^ 1.7i?o in this case, is also reached at similar time {t ^2ms). However, instead of coming to 
rest, a strong retraction of the fluid is seen immediately after reaching and eventually drops to back to zero at 
t = 5.25ms, which corresponds to the moment when the drop completely detaches from the substrate (see Fig. 5b). 

Drop rebound requires that a large portion of surface energy be converted back into kinetic energy during the 
retraction stage. Thus, even though the system is still in the capillary regime, a small amount of viscous dissipation 
or friction can cause the retraction speed to decay, preempting detachment. We therefore expect there to be an 
upper limit for the spreading radius or, equivalcntly, for the impact velocity, above which dissipation becomes large 
enough to eliminate rebound. 

To test this, we carried out a series of experiments with drops of 0.2M HCl washed, unoxidized eCaln, parame- 
terizing the impacts by the Weber number. Fig. 6 shows traces of R{t)/RQ with fixed Rq ^ 0.98mm as function of 
time for different initial Weber numbers {We — 5.6, 19.8, 22.6, 43.2). Since all samples were sufhciently acid-washed, 
the skin effect does not need to be considered. 
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FIG. 6: Weber number dependence of R{t)/Ro for liquid eGaIn washed by 0.2M HCl. The initial drop was fixed at Rq = 
0.98mm. The impact velocity is then chosen to be Vb = 0.45m/s, 0.83ms, 0.90m/s and 1.25m/s. As a result, the corresponding 
Weber numbers are We = 43.2, 22.6, 19.8 and 5.6. 



Clearly, the rebound has disappeared for We — 22.6 and 43.2. Instead, without losing the contact of the substrate, 
a weak retraction follows the maximum radius (t ^ 1.9ms). After the receding velocity vanishes (t ^ 4.7ms), the 
remaining kinetic energy causes the inertial oscillation of the contact line. Finally, the radius stops moving when all 
the kinetic energy is lost from viscous dissipation. On the other hand, at relatively low Weber numbers. We — 5.60 
and 19.8, the detachment between the drop and surface shows up as expected. Therefore, the critical Weber number 
Wee, indicating the upper limit of the rebound regime, is around 20. 

In fact, there should be also a lower Weber number limit to observe rebound, since kinetic energy may not be large 
enough to lift up the drop if the impact velocity is too small. However, in order to explore this regime, the drop 
height has to be so small (<C 1cm) that residual fluctuations in drop velocity as well as drop size, unavoidable with 
our set-up, give rise to large uncertainties. Thus, we were not able to resolve any lower bound from our experiments. 



B. Scaling of the Deformation 



The discussion so far proved qualitatively that the skin can alter the impact behavior in the capillary regime. In 
this section, a more quantitative description of the role of the surface elasticity is used as a one step further toward 
developing a scaling of the maximum spreading factor Pm = Rm/ Ro with impact parameters. For Newtonian fluids, 
Pm is known to scale as We^l^, reflecting a momentum balance between inertial and surface surface tensionl^^l. 
The universality of this scaling has been tested across a wide range of materials, including water, alcohol, viscous 
glycerin and liquid mercury, allowed for large changes in the intrinsic parameters such as density, viscosity and 
surface tensionl^^l . Nevertheless, none of these experiments involved a situation where a portion of surface energy is 
elastically stored in a surface skin. 

In our experiments we controlled the impact by varying velocity Vq, nozzle radius R and oxidation degree. The 
impact velocity was kept under 2m/ s. Therefore, the ratio of viscous to surface energy E^^/E^ — 3CaP^/8 < ClS'^'^l, 
so that viscous dissipation was much less important in this regime (here Ca — nv/'y is the Capillary number). 
Meanwhile, we used nozzles with two radii, R = G.48mm and 4.1mm, to test for size-dependence. 

Fig. 7a plots Pjn vs. We under these conditions. Different colors stand for different surface oxidation levels. The 
solid symbols indicate data taken with the small nozzle while the open symbols correspond to the large nozzle. For 
the eGaIn sample washed with 0.2M HCl before impact to remove any skin (red data points), the data for both 
nozzle sizes collapses onto the conventional We^l'^ scaling (dashed line). This confirms that pure eGaIn follows the 
behavior of Newtonian fluids, for which the scaling of Pm with We should be independent of the drop size. Similar 



9 



results are also observed for drops washed with O.OIM HCl (pink), except that the entire data set shifts downward 
by about 5%. 

However, the spreading of oxidized eGaIn exhibits a significant dependence on the size of the nozzle. For a given 
Weber number We, Pm for the large nozzle (open blue symbols) is seen to fall below that for the small nozzle 
(solid blue symbols), shifting downward by approximately 60%. This is inconsistent with the properties expected of 
Newtonian fluids. We note that for the oxidized eGaIn drops from the 0.48mm nozzle, Rm/Ro remains consistent 
with the conventional scaling Pm ~ We^^^ and significant variance of Pm occurs only for the 4.1mm nozzle. This 
suggests that the skin effect becomes important only when the drop size gets sufficiently large. 
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FIG. 7: (a). Spreading factor Pm vs. Weber number We (a) and effective Weber number We*{h) under different experimental 
conditions. Colors (red, blue, pink, green) indicate experiments performed at different oxidation conditions. Data for two 
nozzle radii, R — 0.48mm (solid circles) and R = 4.1mm (open circles), are shown. The inset of (b) shows the diameter 
l{= 2Ro) of an eGaIn drop. 

Based on the comparisons between the drop mass and Rq, the uncertainty introduced by non-spherical shape is 
less than 8%, which is not large enough to give rise to the scatter of data in Fig. 7 (a). Instead, This drop size 
dependence of Pm can be explained by the rheology results discussed in Sec. III. If the sample is oxidized, the 
skin induces extra in-plane stress (r^) in the skin surface. As a consequence, the effective surface tension includes 
contributions from both the bulk (cr) and the skin (r^) and can be written as (Jeff — cr + Ts = a + G'^^s- As discussed 
in Sec. Ill, G'^ is connected to bulk modulus G' by multiplying with the size scale of the eGaIn sample Zc(which is 
the tool radius for the rheology measurement), G'^ = (G" • ^c)/4. Thus, 

CTe// = cr + = O- + G'^Js 

= a+^{G'-l,hs^a+^{G'-j,)L (5) 

Since there is no relative displacement between the bulk and skin (7 = 7s) and the measured stress of oxidized eGaIn 
directly gives the yield stress Ty, G' -js = G' -j = Ty. Accordingly, equation (5) can be expressed as CTg// = ct + Tj,^c/4. 
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Previous experimentt^l on the surface tension of different sizes eGaIn drops quantitatively proved that the size scale 
Ic of pendant drop is simply the diameter / = 2Ro. As a result, the effective surface tension scales as 



<^eff =(^+^ryRo. (6) 

This model fits the apparent surface tension (cr) from the pendant drop sizes'^' . 

Physically, the size dependence of Uejf comes from the increase of with the drop size. However, both r, and 7^ 
are difficult to directly measure in the experiment.. By relating the surface parameters (Gg,7s) with the measured 
bulk modulus (G") and yield stress (r^) in Eq. (5), we are be able to quantitatively calculate aeff via Eq. (6). 

If we now account for the surface energy stored in the skin by replacing a in We with agff given by Eq. (6) we 
can define an effective Weber number as 

We^='^= (7) 
aeff l + TyRo/2a ^ ' 

The native surface tension of pure liquid eGaIn is cr 4 x lO^mN/m. Pre-washing the sample with 0.2M HCl reduces 
the yield stress to O.lPa (Fig.2). Since the largest drop made in the experiment has the cross diameter Rq ~ 6mm, 
TyRo/{2a) < 0.01 ^ 1 and hence We* ~ We. Eq. (7) also explains why the data for small, oxidized eGaIn drops 
follows the conventional scaling even though the skin is present. However, when the product of I and Ty becomes 
large, We* can differ significantly from We. For instance, with yield stress in air Ty ~ lO^Pa for an oxidized drop 
of size Ro ^ 5mm, TyRo/2a ^ 1 and We* drops below We (open blue symbols in Fig. 7a). To test the validity of 
Eq. (7), we replot from Fig. 7a against the rescaled Weber number We* in panel (b). Within the experimental 
uncertainty, the data collapse nicely onto Pm ^ We*^/^. 

Since the impact velocity was below 2m/s in these experiments, the impact shear stress was below the order of 
lO^Pa. In this regime the loss modulus of the skin is smaller than the elastic modulus by two orders of magnitude as 
mentioned earlier. Hence, viscous dissipation in the skin during impact is negligible and this is why the role of the 
skin can be represented simply by an effective Weber number. 



VI. VISCOUS REGIME 

The surface tension of liquid eGaIn is about one order of magnitude larger than in normal fluids. Therefore, the 
liquid/gas interfacial energy is more important than viscous dissipation in most impact experiments. It is also the 
reason why most previous experimental data for liquid metalsI^^'H^^l'P^l was restricted to the capillary regime. 

In order to observe the role of viscous damping, the impact velocity has to be increased such that a large part 
of initial energy is dissipated through viscous resistance. Experimentally, faster impact was achieved by mounting 
the syringe pump at the top of a long steel rail to increase the drop release height to as large as two meters. This 
extended the impact velocity range up to 6.3m/s. 

In Fig. 8 we compare results for eGaIn drops from the large nozzle (i? = 4.8mm), one air oxidized and another 

pre-washed with 0.2M HCl. The diameter of the drops was (12.3 ± 0.5)mm. Fig. 8a displays the image sequences 

during the impinging of the two drops at impact velocity 3.8m/s. Despite the difference in the initial drop shape, 

the maximum spreading diameter Pm, approximately stays the same. This independence of Pm is also found when 

altering the surface oxidation conditions. In Fig. 8b, Pto(O) indicates the maximum spreading factor of unwashed 

eGaIn drops while Pm{CHCi) represents that of drops washed in an acid bath of specified concentration. Varying 

Chci oyci three orders of magnitude, from O.OOIM to l.OM, the ratio of PmiCHCi) to P„,(0) is seen to fluctuate 

around unity within the experimental uncertainty. This demonstrates that the skin effect is not important in this 

regime and most viscous dissipation is caused in the bulk rather than by the surface. 

n fii 

For Newtonian fluids, the spreading factor scales as Pm ^ 

i?eV5'^°J in the viscous regime. Here the Reynolds 
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FIG. 8: (a). Image sequence of impact process of liquid eGaIn in the viscous regime. In particular, the images on the left 
and right column are taken from oxidized and pure eGaIn respectively, (b). the ratio of PmiCnci) / Pm{0) is plotted against 
the acid concentrations of acid solution used to pre- wash the sample, (c). The classic one-fifth power law can be applied to 
the data of Pm vs. Re in the viscous regime. The Reynolds number is varied by changing the impact velocity and nozzle size 
(4.1mm and 0.48mm). 

number indicates the ratio of inertial to viscous stresses and is defined as 

Re = ^ (8) 
A* 

Given that the skin hardly affects the behavior in the viscous regime, we expect this scaling to survive even for 
oxidized drops. Indeed, by applying a power law fit to the data within more than one decade of Re and different 
oxidation conditions (Fig. 8(c)), we obtain Pm — 0.41i?e^/^. It is worth mentioning that, even though the Reynolds 
number here is around the order of 10^, we do not consider the dissipation caused by turbulence since the time scale 
of impact (^ a couple of milliseconds) is too short to form any inertial flow in the fluid. 



12 



VII. CAPILLARY TO VISCOUS TRANSITION 




FIG. 9: Capillary-to-viscous transition for the impact behavior of eGaIn drops pre-washed in different acid concentrations. 
The dimensionless parameters Pm/Re^^^ and K = We* / Re^^^ are used to collapse all data. 

From the scaling behavior discussed in the last two sections, we can assemble a phase diagram that delineates the 
transition from capillary to viscous regimes. To this end we utilized the renormalized spreading factor Pm/ Re.^^^ and 
the impact number K = We* / Re^/^ X^^"^ Plotted log-log in terms of these variables, in the capillary regime, where 
Pm We*^l^, the data lie on a line of slope 1/4, i.e., Pm/Re^^^ ^ K^^'^. In the viscous regime, where Pm ^ Re^/^^ 
the data lie on a horizontal line, i.e., Pm/Re^^^ = const. 

As seen from Fig. 9, this way of plotting produces excellent collapse of all our data across the full range of acid 
concentrations, drop heights and nozzle diameters. It also highlights the transition from capillary to viscous regime 
at a critical impact parameter value Kc ~ 0.1. The observed value is smaller than in Newtonian fluids, such as 
oil and water, for which ^ l.'^^l This may be attributed to the factors of the surface oxidation or drop shapes 
(Sec. IV), which are not included in the scaling arguments of ordinary fluids. 



VIII. CONCLUSIONS 

In this paper, we experimentally investigated the influence of surface oxidation on the impact dynamics of eGaIn 
liquid metal drops. We showed that the large yield stress in eGaIn exposed to air is associated with an oxide skin 
that can be tuned, as well as eliminated, by immersing the sample into acid solutions of varying concentration (Fig. 
2). In the presence of the skin, the shear modulus G' of liquid eGaIn displays a sample size dependence (Fig. 3a), 
which indicates that G' is not an appropriate intrinsic parameter. Instead, a rescaled surface modulus G'^ — (C ■t)/A 
characterizes the elastic response (Fig. 3b). 

Skin-induced stresses also affect drop impact and spreading. In the capillary regime, where the surface energy 
dominates the dynamics, the oxidation level determines the impact outcome dramatically. At low Weber number, 
unoxidized liquid eGaIn drop can rebound while oxidized eGaIn drops always stay on the surface (Figs. 5 and 6). 
Rebound is enabled by the fact that chemical reduction of the oxide layer can turn the metal surface hydrophobic. 

The skin effect directly affects the geometric deformation during the impact and the subsequent lateral spreading. 
For skin-covered liquid eGaIn in the capillary regime, the spreading factor Pm does not follow the Weber number 
scaling Pm ^ We^^^ observed for ordinary liquids (Fig. 7a). We attribute this to the skin applying an extra line 
tension at the liquidair interface, resulting in an effective surface tension that combines contributions from the 
native material and the oxide skin. We show that this can be captured by an effective Weber number We* = 
We/{1 + Tyl/{4:a)) which collapses all spreading data in the capillary regime when plotted as Pm ~ We*^/'^ (Fig. 
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7b). In particular, from the definition of We* wc sec that the influence of the skin increases with drop diameter I 
and that Pm ^ We*^l^ gives a condition for the effect to become significant. 

Viscous dissipation becomes a dominating factor only at large impact velocities. In this regime, we found no 
difference in spreading factor between pure and oxidized eGaIn (Fig. 8) and the usual Reynolds number scaling 
Pjn ^ Re^/^ is still applicable. Therefore, the resistance to spreading is mostly due to viscous drag in the bulk of the 
material and has very little to do with the skin. Finally, when the renormalized spreading factor Pm/Re^^^ is plotted 
as a function of impact number K = We* / Re^/^ all data obtained under different oxidation conditions, different 
nozzle diameters and different impact velocities, collapses onto a single graph, with a transition between capillary 
and viscous regimes K = We* / Re^/^ ~ 0.1. 

Nevertheless, the importance of the skin-induced surface stress is usually neglected in modeling the impact dynamics 
of molten metal drops'^"''''^^!'!^*' . In these models, the maximum spreading factor Pm was directly obtained from the 
balance between kinetic, native surface energy and viscous dissipation. However, it turned to be difficult to predict 
the spreading radius, especially in the low Weber number regime, without the loss of gener ality[i4].[i5,][i8].Our findings 
suggest the cause of the difficulties is neglecting the energy stored in the oxide skin. Taking its effect into account 
through an appropriately rescaled Weber number We* gives excellent predictability for P^ (Fig. 9). 

Finally, our results may also provide new insights about the impact behavior of general yield stress fluids, including 
dense suspension and polymer gclsl^^l'l^^l . For these classic yield stress fluids, surface properties remain similar to 
ordinary liquids while the kinematic viscosity usually displays unconventional character. Therefore, their unique 
properties should be observed mostly in the viscous regime, which is opposite to the impact of oxidized liquid metals. 

We thank Wendy Zhang and Michelle DriscoU for many useful discussions, and Qiti Guo for the help of preparing 
the experimental samples. This work was supported by the National Science Foundation (NSF) MRSEC program 
under DMR-0820054. 
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